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We consider NLO chargino production and decays at the ILC. For this, we 
present an NLO extension of the Monte Carlo Event Generator WHIZARD in- 
cluding the NLO production. For photonic corrections, we use both a fixed 
order and a resummation approach. The latter method evades the problem 
of negative event weights and automatically includes leading higher order 
corrections. We present results for cross sections and event generation for 
both methods. As a first step towards a full NLO Monte Carlo, we consider 
a LO implementation of the chargino production and subsequent leptonic 
decay and investigate the precision of the sneutrino mass determination 
by means of lepton energy distributions in chargino decays. The SM and 
SUSY backgrounds are included in our study using full matrix elements as 
well as smearing effects from ISR and beamstrahlung. Without using en- 
ergy distribution fits, the sneutrino mass can be determined with an error 
in the percent regime. 
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1. Introduction 

In many GUT models, the masses of charginos tend to be near the lower 
edge of the superpartner spectrum, and can be pair-produced at a first- 
phase ILC with cm. energy of 500 GeV. The precise measurement of their 
parameters (masses, mixings, and couplings) is a key for uncovering the 
fundamental properties of the Minimal Supersymmetric Standard Model 
(MSSM) pp. Regarding the experimental precision which is in the percent 
regime at the ILC, off-shell kinematics for the signal process, and the re- 
ducible and irreducible backgrounds [2] need to be included as well as NLO 
corrections for chargino production. Here we present the inclusion of the 
latter [31 Q]. For decay modes, we focus on the leptonic decay with electron 
and muon in the final state. If sneutrinos decay invisibly into the LSP and 
a neutrino, this channel provides tools to determine sneutrino masses. Such 
decays, common in many scenarios within the MSSM, preclude threshold 
scans since sneutrinos cannot be reconstructed directly. The only possibil- 
ity to access the sneutrino mass in such a case is to select a cascade decay 
where the precise determination of kinematic distributions gives access to 
the sneutrino mass. Although this idea has already been exploited [5], a 
thorough study of how precise the mass determination for the sneutrino in 
the environment of an ILC can be has as yet not been made. We study the 
pollution effects of all reducible and irreducible SM and SUSY backgrounds 
on the visibility of the signal as well as the precision of the sneutrino and 
the chargino mass measurements. We restrict ourselves to areas in SUSY 
parameter space where charginos are within reach of a 500 GeV ILC. 

2. Chargino production at LO and NLO 

Fixed order approach 
The total fixed-order NLO cross section is given by 

O"tot0, m 2 e ) = O-Born(s) + o- v+s (s, A^.^e) + ^-^(s, AE y , ml), 

where s is the cm. energy, m e the electron mass, and A E~ the soft photon 
energy cut dividing the photon phase space. The 'virtual' contribution 
a v is the interference of the one-loop corrections [B] with the Born term. 
The collinear and infrared singularities are regulated by m e and the photon 
mass A, respectively. The dependence on A is eliminated by adding the 
soft real photon contribution a s = / so ft CBom('S) with a universal soft factor 
/soft(-^x 1 ) [7J. We break the 'hard' contribution cr 2 _»3(s, A-Ey, rn%), i.e. the 
real-radiation process e~e + — ► XiXjl-i a collinear and a non-collinear 
part, separated at a photon acollinearity angle A# 7 relative to the incoming 
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electron or positron. The collinear part is approximated by convoluting the 

2 

Born cross section with a structure function fh(%] A# 7 , ^-) [8]: 

rEmax,A9 

o- h , c (AE,A9,s) = / dx i dr 2 f h (x l )\M h \ 2 {x l ,s), (1) 

JAE.O 

where Xi denotes the momentum fraction of the respective incoming beam 
after photon radiation and (W2 the two particle final state phase space. The 
non-collinear part is generated explicitly using exact three particle final state 
kinematics. 

The total fixed order cross section is implemented in the multi-purpose 
event generator WHIZARD [9l [10] using a 'user-defined' structure function 
and an effective matrix element 

\M cS \ 2 = (l + / soft (A£ 7 , A))|.M B orn| 2 + 2Re(M Bm nM; irt (X)) 

which contains the Born part, the soft-photon factor and the Born-1 loop 
interference term. In the soft-photon region this approach runs into the 
problem of negative event weights [11]: for some values of 9, the 2^2 
part of the NLO-corrected squared matrix element is positive definite by 
itself only if AE^ is sufficiently large. To obtain unweighted event samples, 
an ad-hoc approach is to simply drop events with negative weights before 
proceeding further. 

Resummation approach 

Negative event weights can be avoided by resumming higher-order initial 
radiation using an exponentiated structure function /isr [121 S3]- I n order to 
avoid double-counting in the combination of the ISR-resummed LO result 
with the additional NLO contributions [BJ, we have subtracted from the 
effective squared matrix element the soft and virtual photonic contributions 
that have already been accounted for in o~ s +v This defines 

\Mlf\ 2 = \M eS \ 2 - 2/ soft , IS R |.MBorn| 2 (2) 

which is positive even for low AEy cuts for all values of 9. Convoluting 
this with the resummed ISR structure function for each incoming beam, we 
obtain a modified 2 — ► 2 part of the total cross section which contains all 
NLO contributions and in addition includes higher order soft and collinear 
photonic corrections to the Born/one-loop interference. This differs from 
the standard treatment in the literature (cf. [6j) where higher order photon 
contributions are combined with the Born term only ("Born+"). 
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Fig. 1. Total and relative cross section as a function of yfs. Left figure: Born (red, 
"LO"), fixed order (blue, "NLO") and fully resummed (green, "NLO") total cross 
section, right figure: relative fixed order (blue, solid) and fully resummed (green, 
dashed) higher corrections with respect to Born result 



The complete result also contains the hard non-collinear 2 
convoluted with the ISR structure function: 



Or. 



A(E,i 



dx i dT 2 f^\x 1 )f^ R '(x 2 )\M^ 



4 res 1 2 



A(E,6) 



dx i dT 3 f^\x 1 )f^ R \x 2 )\M 2 ^ 3 \ 2 



3 part 



(3) 



The resummation approach eliminates the problem of negative weights such 
that unweighting of generated events and realistic simulation at NLO are 
now possible in all regions of phase-space. 



3. NLO Chargino Production: Results 

Total cross section and relative corrections 

Figure [T] shows the cm. energy dependence of the total LO and NLO 
cross section for chargino production for the mSugra point SPSla' [T] and 
the relative corrections with respect to the Born result. The corrections are 
mostly in the percent regime and can reach 20% in the threshold region. 



Cutoff dependencies 

Figure [2] compares the AE^ dependence of the numerical results from a 
semianalytic fixed-order calculation with the Monte-Carlo integration in the 
fixed-order and in the resummation schemes. The fixed-order Monte-Carlo 
result agrees with the semianalytic result as long as the cutoff is greater 
than a few GeV but departs from it for smaller cutoff values because here, 
in some parts of phase space, |-M e ff| 2 < is set to zero. The semianalytic 
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Fig. 2. Total cross section dependence on A£" 7 : 'sa' (red, dotted) = fixed-order 
semianalytic result; 'fix' (green, dashed) = fixed-order Monte-Carlo result; 'res' 
(blue, long-dashed) = ISR-resummed Monte-Carlo result; (magenta, dash-dotted) 
= same but resummation applied only to the 2 — ► 2 part. A0 7 = 1°. LO: Born 
cross section. 



fixed-order result is not exactly cutoff- independent, but exhibits a slight rise 
of the calculated cross section with increasing cutoff due to the breakdown 
of the soft photon approximation. For AE^ = 1 GeV (10 GeV) the shift 
is about 2 %o (5 %o) of the total cross section. The fully resummed result 
shows an increase of about 5 %o of the total cross section with respect to the 
fixed-order result which stays roughly constant until AEy > 10 GeV. This 
is due to higher-order photon radiation. 

For the dependence on the collinear cutoff A# 7 , the main higher-order 
effect is associated with photon emission angles below 0.1°. For A# 7 > 10°, 
the collinear approximation breaks down. 



Event distributions 

In Fig. [3] we show the binned distribution of the chargino production 
angle obtained using a sample of unweighted events. It demonstrates that 
NLO corrections to the angular distribution are statistically significant and 
cannot be accounted for by a constant K factor. 



4. LO production and leptonic decays 

Signal and (MS)SM backgrounds 

In the second part of our work, we investigate the LO chargino produc- 
tion and subsequent leptonic decay modes. To avoid large SM backgrounds 
arising in the production of same flavor and opposite sign lepton pairs, we 
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Fig. 3. Polar scattering angle distribution for an integrated luminosity of 1 ab _1 
at y/s = 1 TeV. Left: total number of events per bin; right: difference w.r.t. the 
Born distribution. LO (black, dotted) = Born cross section without ISR; fix (green, 
dashed) = fixed-order approach; res (blue, full) = resummation approach. 



consider opposite flavor opposite sign lepton pairs in the final state: 

e + e" — > xtXi — ► XiX?e~ /U + v^v^. (4) 

This channel is especially interesting when the sneutrinos decay invisibly 
and the sneutrino pair-production channel is therefore experimentally inac- 
cessible. In chargino pair-production and subsequent decays, however, the 
sneutrino and chargino masses can be determined from the edges of the 
lepton energy distributions [5] . The experimental signature of the signal ^ 
is 

e+e~ — ► e"/i + + E miss . 

Therefore, all background processes where the missing energy results from 
the emission of invisible particles (z/s, x°' s ) need to be considered. In 
addition, we have to take into account processes where additional particles 
are emitted at very small angles and vanish in the beampipe. We therefore 
investigate the SM and MSSM backgrounds as listed in Table 1. 

The values of the total cross sections for the SUSY parameter point 
SPSla' including beamstrahlung and initial state radiation before the ap- 
plication of cutgj are given in Table [2j The most dominant background 
process is photon induced t + t~ production; its cross section exceeds the 
magnitude of the signal cross section by a factor 10 4 . Similarly, SM back- 
ground processes such as direct W and r (pair) production are significantly 
larger than the signal, while SUSY backgrounds are of similar size. 



For numerical reasons, we always include a collinear cut of 5° for the outgoing e 
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ID 


final state 


most dominant process 




e + e~e~ v e v T v r 


7-induced r pair production (SM) 


WW 




WW production (SM) 


T 




r pair production (SM) 


tW 


e~ jj, + v^ v e v T v T v T u T 


r from WW production (SM) 


jW 


e _ e + e~/x + f M v e v T v T 


7-induced WW production (SM) 


f 




f pair production (MSSM) 




e'ti+Vp V e V T V T V T U T XlXl 


tv t from x decays (MSSM) 



Table 1. SM and MSSM background processes leading to e /j + + E 1 



ID 


before cuts 


after cuts 


signal 


3.940 (8) 


1.905 (4) 




25495 (4) 


0.072(1) 


WW 


152.42 (41) 


0.794(2) 


T 


34.8(18) 


0.024(1) 


tW 


2.978 (9) 


0.185(1) 


jW 


2.192(12) 


0.154(1) 


f 


4.107(7) 


1.146(2) 


f v T 


2.74(9) 


0.72(2) 



Table 2. Signal and background total cross sections before and after the appli- 
cation of background suppression cuts, for SPSla' and y/s — 500 GeV. ISR and 
beamstrahlung included. All results are given in fb, together with an integration 
error. 



p±(e,fi) > 2 GeV, p±(e) + P± (jj.) > 4 GeV, 
lGeV < E(e,fj) < 40GeV, -160° < < 160°, 
15° < 6{e) < 155°, 25° < 9{ji) < 165° 

Table 3. Cuts applied for background suppression. A(j) is the azimuthal separation 
angle of the lepton pair. 

5. LO production and decay: Results including cuts 

In order to suppress the SM and MSSM background, we apply the set 
of cuts given in Table [3l The magnitudes for the total cross sections after 
the application of these cuts are presented in Table [2 The signal has been 
reduced by roughly a factor 2, while the dominant background from photon 
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Fig. 4. Energy distribution of electron (left) and muon (right) after cuts given 
in Table EH Grey: SM background; (photon-induced) t + t~ and W + W~ pair 
production. Red: SUSY background; ff and ttv t v t production. Blue: SUSY 
signal; xfxi production and successive leptonic decays. 



induced tau pair-production is now suppressed by 10 6 . In total, we obtain 
a signal/background ratio of 0.62, which, for an integrated luminosity of 
lab -1 , leads to a 20a discovery. Figure [4] shows the energy distribution 
of the leptons after cuts have been applied. While the SM processes lead 
to a flat background distribution which can be easily subtracted, SUSY 
background processes are more challenging, as they result in kinematic dis- 
tributions similar to the signal. 

Chargino and sneutrino mass determination 

In the case of quasi mass-degenerate sneutrinos and the lightest charginos, 
the leptonic decay mode is dominated by intermediate on-shell sneutrinos: 

x± -» l±v x - /±^x°- (5) 

In these scenarios, sneutrino masses can only be determined from the edges 
in the lepton energy distributions in chargino decays ©. From on-shell 
relations, we obtain [5] 

i— V^tnin -E'max /-, 2 (-Emm ~\~ ^max) , n \ 

m x± = — Tf — ' 1719 = "^V 1 7f ' (6) 

-C'min ~r "max y V s 

where -Emm/max are the minimum/maximum lepton energy. Taking naive 
readoff values for the energies, we have (cf. Fig. 0|) 



E min = 4.5 ± 1.0 GeV, E max = 24.5 ± 2.0 GeV 
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and, using ([6]), 

m^± = 181 ± 15GeV (183.67), mo = 170 ± 14GeV (173.52), 

where the values in brackets are the nominal (input) values for SPSla'. 
Although the central values are in good agreement, the large errors here 
call for a refined treatment. Alternatively, we take the chargino mass from 
threshold scans; assuming A t h r f^x = ^ GeV, we then obtain 

mc = 172.7 ± 1.3 GeV, 

where the error is now in the percent regime. A fitting routine for error 
reduction down to a permille range using edge distributions including all 
backgrounds is in the line of future work. 

6. Conclusions 

We have implemented NLO corrections into the event generator WHIZARD 
for chargino pair-production at the ILC with several approaches for the in- 
clusion of photon radiation. A careful analysis of the dependence on the 
cuts AEy, A9y reveals uncertainties related to higher-order radiation and 
breakdown of the soft or collinear approximations. Careful choice of the 
resummation method and cutoffs will be critical for a truly precise analysis 
of real ILC data. The version of the program resumming photons allows to 
get rid of negative event weights, accounts for all yet known higher-order ef- 
fects, allows for cutoffs small enough that soft- and collinear-approximation 
artifacts are negligible, and explicitly generates photons where they can be 
resolved experimentally. 

Additionally, we have investigated the complete production and decay 
process at LO including several SM and MSSM backgrounds. For this, we 
used a full matrix element simulation including initial state radiation and 
beamstrahlung. Efficient cuts reduce the dominant background by a factor 
10 6 allowing a clear isolation of the signal. The leptonic decay mode allows 
for mass determination of sneutrinos from the edges of the lepton energy 
distributions. Taking the chargino mass from threshold scans, the error 
of the sneutrino mass determination can be reduced to the percent level. 
The improvement of the mass determination precision using refined fitting 
routines and the combination of NLO production and decay |14t [T5] are in 
the line of future work. 
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